Phonon bottleneck in the low-excitation limit 
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The phonon-bottleneck problem in the relaxation of two-level systems (spins) via direct phonon 
processes is considered numerically in the weak-excitation limit where the Schrodinger equation 
for the spin- phonon system simplifies. The solution for the relaxing spin excitation pit), emitted 
phonons rik(t), etc. is obtained in terms of the exact many-body eigenstates. In the absence 
of phonon damping F p h and inhomogeneous broadening, p(t) approaches the bottleneck plateau 
Poo > with strongly damped oscillations, the frequency being related to the spin-phonon splitting 
A at the avoided crossing. For any r p h > one has p(t) — » but in the case of strong bottleneck 
the spin relaxation rate is much smaller than F p h and p(t) is nonexponential. Inhomogeneous 
broadening exceeding A partially alleviates the bottleneck and removes oscillations of p(t). The line 
width of emitted phonons, as well as A, increase with the strength of the bottleneck, i.e., with the 
concentration of spins. 

PACS numbers: 31.70.Hq, 63.20.-e, 67.57.Lm 



I. INTRODUCTION 

Spin-lattice relaxation is an old an much studied prob- 
lem that currently recieves a resurge of attention because 
of its vital importance in quantum information process- 
ing (see, e.g., Ref. 0] and references therein). Theoretical 
description of the spin-lattice relaxation as a single-spin 
process is in many cases insufficient because of the col- 
lective effects of incoherent and coherent nature, such as 
the phonon bottleneck 2 and superradiance,- respectively. 

The problem of phonon bottleneck (PB) in relaxation 
of two-level systems (henceforth spins) via direct phonon 
emission/absorption processes, first recognized by Van 
Vlcck 2 in 1941, remains unsolved until now. In two 
words, if the emitted phonons have nowhere to go, they 
are absorbed by spins again and thus the spins cannot 
relax efficiently. However transparent this picture might 
appear, is not easy to propose a theoretical description 
of the effect based on the first principles. 

Published theories of the PB2£&1& use ad hoc rate 
equations for populations of spins and resonant phonons, 
considering the latter as a single dynamical variable. 
This is certainly an oversimplification, because the emit- 
ted phonons, having frequencies Wk, form a group with 
a bell-like line shape with some width, centered around 
the spin transition frequency luq. For a single spin em- 
bedded into an infinite elastic matrix (as well as for a 
decaying atomic state in a free space) this line shape is 
Lorentzian with the width r/2, where T is the single- 
spin decay rate following from the Fermi golden rule^^ 
However, in the case of many spins with a concentration 
sufficient to create a bottleneck, the line shape and line 
width of emitted phonons are unknown and should follow 
from the solution of the problem. 

Van Vleck came to the idea of the phonon bottle- 
neck comparing the rate of energy transfer from spins 
to phonons (obtained using experimental data) with the 
phonon relaxation rate T p h due to different mechanisms 
and he found the latter to be typically too small to keep 
the phonon subsystem at equilibrium. However, the pri- 



mary role in the PB problem belongs to another parame- 
ter that is not related to the phonon relaxation rate. This 
parameter is of a statistical origin and is defined as the 
ratio of the number of spins to the number of phonon 
modes within the single-spin line width T. 11 If this so- 
called bottleneck parameter B is vanishingly small, the 
spin excitation goes over into the phonon subsystem and 
never returns. In this case spins completely relax even 
without any phonon damping. However, for nonzero B 
and T p h = 0, the spin relaxation ends in the so-called 
bottleneck plateau that corresponds to a quasiequilib- 
rium between spins and resonant phonons but not to the 
complete equilibrium. Further relaxation to the com- 
plete equilubrium can be achieved only if T p h is taken 
into account. It should be stressed that the effective re- 
laxation rate of the spins in this case is not r p h. It is 
much smaller and can be estimated as T p h multiplied by 
the small fraction of phonon modes in the total number 
of modes (phonons + spins) involved in the process. 

Although in many practical situations the number of 
resonant phonon modes is determined by the inhomo- 
geneous broadening of spin levels, the pure spin-phonon 
model without inhomogeneous broadening has a funda- 
mental importance. It was shownli that this model can- 
not be described kinetically (i.e., in terms of spin and 
phonon populations only) because of long-memory ef- 
fects. In Ref. Oj] memory effects have been taken into 
account within a minimal approximation, adding a new 
variable that can be interpreted as spin-phonon correla- 
tor. Analytical and numerical solutions of the resulting 
reversible dynamical equations show that the spin ex- 
citation approaches the botleneck plateau with damped 
oscillations. Inclusion of an ad hoc phonon damping r p h 
into the bottleneck equations allows to describe the sec- 
ond stage of the relaxation towards the complete equilib- 
rium. 

Still, the solution of the PB problem in Ref. [ll| is not 
completely satisfactory since it cannot produce a well- 
behaved line shape of emitted phonons. This indicates 
that additional nondiagonal correlators should be taken 
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into account that will make the description more com- 
plicated. Another important factor that should be taken 
into account is the inhomogeneous broadening of the spin 
levels. 

This paper presents the exact numerical solution of 
the phonon bottleneck problem based on the Schrodinger 
equation for the spin-phonon system, with and without 
the ad hoc phonon damping and inhomogeneous spin 
broadening. The full Schrodinger equation for a many- 
body system is, of course, intractable by direct methods 
because of too many variables. However, the low-excited 
states of the system can be described by a single excita- 
tion that is hopping between spins and phonon modes. 
In this case the Hilbert space of the problem is severely 
truncated and one has to work with matrices the size 
of which is just the total number of spins and phonon 
modes under consideration. This is the case considered 
here and the solution of the PB problem is obtained by 
matrix algebra using Wolfram Mathematica. The results 
of the calculations show that for the pure model the spin 
excitation approaches the bottleneck plateau with oscil- 
lations, however, less revealed than in Ref. [nj. The line 
width of emitted phonons broadens with the bottleneck 
parameter B. 

The structure of the rest of the paper is the following. 
Sec.|TT]sets up the Hamiltonian of the spin-phonon system 
and the Schrodinger equation (SE) in the case of a single 
excitation. The spin excitation and the initial conditions 
for problem are defined here. Sec. IIIII presents the known 
results for the relaxation of a single spin and for the en- 
ergy distribution of emitted phonons, used later for the 
reference. Sec. IIVI introduces the bottleneck parameter 
B for systems of many spins from statistical arguments, 
both with and without inhomogeneous spin broadening. 
Sec. [V] is the central section of the paper introducing the 
matrix formalism for the single-excitation spin-phonon 
problem. Here the expressions for the spin excitation 
p(t), its asymptotic value p^ (the bottleneck plateau), 
and the asymptotic populations of the emitted phonons 
are obtained in terms of eigenvectors and eigenvalues of 
the dynamical matrix of the system. General analysis of 
the eigenstates of the dynamical matrix is done in Sec. 
IVII It is shown here that the number of the phonon 
modes "on speaking terms" with spins for the pure prob- 
lem increases with B thus changing the statistical balance 
between spins and phonon modes. In this section the for- 
mulas describing the spectrum of the split spin-phonon 
modes and the hybridization of different phonon modes 
with each other (i.e., the scattering of phonons on spins) 
are obtained in terms of the eigenstates of the dynamical 
matrix. Sec. lVIII presents the results of numerical calcula- 
tions for the pure model, including the split spin-phonon 
modes with the gap, time evolution of the spin excita- 
tion, bottleneck plateau, the post-plateau relaxation due 
to the phonon damping, and the energy distribution of 
the emitted phonons in the pure model. In Sec. IVIIil thc 
effects of the inhomogeneous spin broadening are consid- 
ered. The latter is shown to wash out oscillations of the 
spin excitation pit). Sec. IIXI contains the implementa- 
tion of the general results to the spin relaxation between 



adjacent spin levels in molecular magnets. In Sec. [X] fur- 
ther problems of collective spin-phonons relaxation are 
discussed. 



II. THE HAMILTONIAN AND SCHRODINGER 
EQUATION 

Consider a spin-phonon Hamiltonian for N$ two-level 
systems (spins) put at positions within an elastic body 
of N cells 



H = H () 



V, 



where 



(1) 



(2) 



describes spins and harmonic phonons, er being the Pauli 
matrix. The spin transition frequencies iOi can differ from 
site to site. One can represent them in the form 



LUQi = U! + 5u} 0l 



(3) 



where Scooi <C is the inhomogeneous broadening. In 
the absense of the latter we will simply use loq as the 
spin transition frequency. Neglecting the processes that 
do not conserve the energy (that can be done in cases of 
practical significance where V can be treated as a per- 
turbation) , one can write V in the rotating- wave approx- 
imation (RWA) as 



7^ E E K^ 01 4 + A tk XlW) , (4) 



V 



where = V^e - r< . In the numerical work below Vk 
will be replaced by a constant, Vk =>■ V. The opera- 
tor A 10 = er_ brings the spin from the ground state 
|T> = |0) = Q to the excited state \[) = |1) = (°) while 
A 01 = <r + does the opposite. Note that the state with 
no phonons and all spins in the ground state is the true 
ground state of the Hamiltonian above. 

Below we will consider the low-excited states of the 
spin-phonon system that can be described by a superpo- 
sition of the vacuum state of the system |0) (no phonons 
and all spins in the ground state) and the states with one 
excitation that is hopping between the spins and phonon 
modes. The wave function of these states has the form 

*= ^o+E c ^ 10 +E c ^ |o>, (5) 



where the coefficients satisfy the system of equations 
dc 



dt 



ILOqCq 



(6) 



and 



dci . . i \ . 

dt ^/N 4^ 

k 

dt 



-i (w k - u a ) c k + -j= A* k 



ck, (7) 
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up to an irrelevant global phase factor. One can see that 
the ground-state coefficient Co is decoupled from the other 
coefficients since the RWA Hamiltonian conserves the ex- 
citation number 



p(t) + nk(f) = const, 



where 



(8) 



(9) 



are the excitation number of the spin subsystem and pop- 
ulations of the phonon modes. One can also define trans- 
verse spin polarization components by 



(*+> = <* 01 ) 

(O = {*+)*, 



c E Ci 



(a x ) = Re (<t + ) , 



(10) 



etc. The main part of the time dependence of (er+) is 
e -Kj t_ -jjjg absolute value of the transverse spin compo- 
nent 

He/wV^HWI (ii) 

does not have this oscillating factor. In many practical 
situations (erj_) decays with time due to the inhomoge- 
neous broadening. In the absence of the latter, the only 
source of the decoherence is interaction with phonons. 

Our task is to find the time evolution p(t) and (cr±) t 
starting from a particular initial state. In this work we 
restrict ourselves to the initial states with no phonons, 
Ck(0) = 0. The simplest initial condition in this case is 
one spin at site io excited and all other spins in their 
ground states: 



Ci (0) = l, c M4o (0) = 0. 



(12) 



Another kind of the initial spin state is the state with the 
excitation equidistributed over all spins: 



c,-(0) - e^/y/Ng. (13) 
The initial spin state with random phases 

(e«**-W) = 5 ij , (14) 

is called incoherent. If cj) i are constant or they periodi- 
cally change in space with some wave vector qo , the initial 
state is coherent. One can consider other kinds of spin 
initial conditions, say, excitation distributed over spins 
in some compact region of space. 



III. NON-BOTTLENECKED SPIN-LATTICE 
RELAXATION 

The results of this section can be found in the 
literature^iiLIi still a concise description of the non- 
bottlenccked spin-phonon dynamics is presented for the 
sake of consistency and future reference. 



A. Relaxation of a single spin 

Suppose there is a single spin, a = c, in the initially 
excited state. With a proper choice of the origin of the co- 
ordinate system one has = Vk. Using the Schrodinger 
equation (J7J), one can integrate the equations for the 
phonon modes c k : 



Ck{t) 



VNJt 

= Ilk, / dre-^'^cU - r) (15) 



and insert the result into the equation for the spin c: 



dc 



1 



t-to 



-di = ~ n E |Fk|2 ' dre-^-^c{t - r). (16) 

k 







In this integro-differential equation, c(t — r) is a slow 
function of time, whereas the memory function /(t) = 



(l/JV)EkWe- 



-i(<*>k— «o)i 



is sharply peaked at r = 0. 



Thus one can replace c(t — r) => c(t) , after which integra- 
tion over r and keeping only real contribution responsible 
for the relaxation yields the equation 



dc 
~dt 



-c, 



where 



r = f El^l 2 ^K-^o) 

k 



(17) 



(18) 



is the single-spin decay rate. For Vk = V independently 
of the direction of k, T can be written as 



r 



27r\V\ 2 p ph M, 



where 



(19) 



(20) 



is the phonon density of states normalized by one. The 
accuracy of the above short-memory approximation is 
justified by T < wo- The ^-function in Eq. (fig)) implies 
that the spin is relaxing to a large number of phonon 
modes so that summation is replaced by integration, 



N 4-^ 



d d k 



(21) 



where d is spatial dimension. In small bodies with es- 
sentially discrete phonon modes Eq. (f2Tj) is invalid. The 
solution of Eq. (TT7J) with t = is c(f) = e~( r / 2 )* that 
leads to well-known decay rule for the spin excitation 



Pit) 



-Ft 



(22) 



The rate of transverse spin relaxation according to Eqs. 
(jTOf and (HU is T/2. 

A similar elimination of phonons can be performed in 
the case of many spins in the absence of the PB^ The 
resulting equations describe collective spin-phonon relax- 
ation, including superradiance^ 
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B. Distribution of emitted phonons 

After the time dependence c(t) has been found, one 
can return to Eq. (Tj"5)) and calculate c q (£). With t = 
the result is 



iV£ e -»(^-"o)< _ e -rt/2 
Ck(t) - k 



/]V -i (w k - w ) + r/2 ' 
This leads to the distribution of emitted phonons 



(23) 



n k (t) = \c k (t)\ A = 



2 |F k | 2 1 - 2e- rt / 2 cos [{u k - lu ) t] + e~ rt 



N 



(wk-wo) +r 2 /4 



(24) 



that asymptotically becomes the Lorentzian function 



rik = 



iV( Wk - Wo ) 2 + r 2 /4 



(25) 



With the help of Eqs. (J22J) and (HU one can check that 
the total excitation is conserved, in accordance with Eq. 
©. For For Vk = V Eq. (|2"5|) can be rewritten as 



n k 



r/2 



irN Pph {u ) ( Wk - Wo ) 2 + r 2 /4 ' 



(26) 



IV. THE PHONON BOTTLENECK 

Let us now turn to systems with a macroscopic number 
of spins Ns- At least in the case of diluted spins, the 
relaxation is controlled by the bottleneck parameter B 
that can be defined as the ratio of the number of spins N$ 
to the number of phonon modes Nr within the natural 
spin line width T of Eq. (fT8)) . 



A i 



^ N P P h ( w o) r, 



where p ph (uj) is given by Eq. (f20|) . That is^ 

N s n s 



B = — 



N r nN Pph (uq)T irp ph (u ) T ' 



(27) 



(28) 



ris = N$/N being the number of spins per unit cell. 
The definitions above pertain to a single phonon branch 
and extension to several phonon branches is obvious. For 
B < 1 (see below) Nr is the estimation of the number of 
phonon modes that can exchange excitation with spins. 

In the case of a single spin, Ns = 1, in a macroscopic 
(N — > oo) matrix the parameter B is vanishingly small. 
The excitation, initially localized at the spin, spreads 
with time over a large number Nr of resonant phonon 
modes, so that the spin relaxes completely according 
to Eq. (j2"2"|) . In simulations, the macroscopic limit is 
achieved if the average distance between the neighbor- 
ing phonon modes becomes smaller than the natural line 
width r 



1 



Np(u ) 



<r. 



(29) 



If the sample is so small that the spin can exchange ex- 
citation with only a few phonon modes, Nr ~ 1, it does 
not relax completely. In this case one has B ~ 1, the 
so-called phonon bottleneck situation. The simplest re- 
alization of the bottleneck is a system of two resonant 
states in which the excitation oscillates between the two 
states in time. If Ns is macroscopic but still Ns <C Nr 
and thus B<1, the initial spin excitation is transferred 
irreversibly into the phonon subsystem, as is clear from 
statistical arguments. There is no bottleneck for B <C 1, 
and the spin relaxation is still described by Eq. (j2"2"]) . 

In the case of a finite concentration of spins ns, the 
parameter B can easily become large. In this case only a 
small fraction of the excitation migrates into the phonon 
subsystem and the spins practically cannot relax as the 
emitted phonons are being absorbed by spins again. For 
B 3> 1, as we will immediately see, the number of 
the phonon modes "on speaking terms" with spins is 
not Nr but much greater. The latter can be obtained 
if one considers the spin-phonon hybridization. Insert- 
ing the spin Fourier components 6k = Cie lk ri into 
the Schrodinger equation ([7]) and neglecting the cou- 
pling of modes with different values of k, i.e., using 
J2i ^e i(k ~ q )' ri =*> N S S (k - q) , one can reduce Eq. © 
to a 2x2 matrix problem that can be easily diagonal- 
ized. The eigenstates of this problem are hybridized spin- 
phonon modes with frequencie a 13 ' 14 



(cj k - w ) + A 2 



where 



A k = 2V^|^k| 



(30) 



(31) 



is the spin-phonon splitting. These modes have the form 
of a straight line with a slope corresponding to phonons 
and a horizontal line corresponding to spins. The lines 
have an avoided level crossing at w k = loo, split by Ak 
that depends on the spin concentration. For Vk = V 
independently of the direction of k, one can eliminate 
p ph (cj ) from Eq. (J2S]) with the help of Eqs. (J201) and 
and obtain the relation 



B 



2ns\yV 

r 2 



2r 2 ' 



(32) 



For A > r the number of phonon modes strongly cou- 
pled to spins can be estimated as 



N A = irNp h (w ) 



(33) 



If ns is so small that A falls below the natural spin line 
width r, one cannot speak of the hybridyzed spin-phonon 
modes. From Eq. (f32|) one obtains 



N A = NrVB/2. 



(34) 



The number of resonant phonons Af res that exchange ex- 
citation with the spins can be estimated in the whole 
range of B as 



iVre 



N r , B<1 
N A , B>1. 



(35) 
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For any B > 0, the spin excitation p(t) does not relax 
to zero but reaches a plateau at some p^ that from the 
statistical equidistribution argument can be estimated as 



Po 



N s 



N s + N Te 



(36) 



with iV r es defined above. We will see that this formula 
works well both for B < 1 and B > 1, in the absence of 
the inhomogeneous broadening that will be considered in 
Sec. I Villi The asymptotes of the above expression are 



B, B < 1 

l-l/VW, S> 1, 



(37) 



in accordance with the numerical results of Sec. I VIII 

To describe the complete spin relaxation after the bot- 
tleneck plateau, one has to include the phonon relaxation 
processes, the easiest way being ascribing an empirical 
relaxation rate T p h to the phonons. 



but they are very slow. An important advantage of this 
method is that it allows to obtain formulas for the asymp- 
totic t — > oo state of the system in terms of matrices. 

Below, the method based on the eigenvectors of <fr will 
be used. In the general case when the ad hoc phonon 
damping r p h is added, 3? is non-Hermitean, and one has 
to distinguish between right and left eigenvectors. The 
dynamical matrix $ has Ns + N right eigenvectors R M 
that satisfy 



(39) 



In the eigenvalues, the imaginary parts originate from 
r p h- The size of the matrix <fr can be reduced if the 
phonon modes far from the resonance are dropped. The 
solution of Eq. (|38|) can be expanded over the complete 



orthonormal set of R M as follows 



C(t) = ^R A1 e-^ +r ^ t L, I • C(0), 



(40) 



V. DYNAMIC MATRIX AND THE TIME 
EVOLUTION OF THE SYSTEM 

For the numerical solution of Eq. ([7]) it is convenient to 
introduce the state vector C = ({ci},{ck}) and rewrite 
Eq. in the form 



dC 
~~dl 



— = C, 



(38) 



where 3> is the dynamical matrix of the spin-phonon sys- 
tem. In the absence of phonon damping, <fr is Hermitean. 
Since the number of discreet phonon modes is N, 4> is a 
(Ns + N) x (Ns + N) matrix. There are three methods 
of numerical solution of this equation that can be imple- 
mented in the Wolfram Mathematica. 

The first method is the direct numerical solution using 
one of the ordinary-differential equations (ODE) solvers. 
This method is fast, can be made much faster if Math- 
ematica is replaced by one of programming languages, 
it does not require high accuracy, but it does not allow 
to analytically average over the random phases in Eq. 
(j!3[) . This averaging can only be done if one runs the 
calculation many times with different realizations of ini- 
tial conditions. 

The second method is based upon numerical calcula- 
tion of the matrix exponentials in the solution of Eq. (j3"5|) 
C(t) = e~ l#t C(0). This method is slower than the direct 
ODE solution, it also does not require high accuracy, and 
here one can average over the initial conditions analyti- 
cally. 

The third method uses the expansion of the solution 
C(t) over eigenvectors of <&. This method allows ana- 
lytical averaging over initial conditions, it is faster than 
the method using matrix exponentials, if formulated in 
a fully vectorized form. However, this method requires 
high precision and for large matrices it runs on 64-bit 
machines only. Note that arbitrary-precision computa- 
tions on 32-bit machines with Mathematica are possible 



where L M are left eigenvectors of <& that satisfy L p • R„ = 
<5^ t „. Note that, in general, R^ t and are not complex 
conjugate. The vectorized form of Eq. (|4"0|) is 



C(t) = E • W(i) • E" 1 • C(0), 



(41) 



where E is the right-eigenvector matrix composed of 
all eigenvectors R M standing vertically, E _1 is left- 
eigenvector matrix, composed of all left eigenvectors lying 
horizontally, and W(t) is the diagonal matrix with the 
elements e - (in " +r f } *. In fact, E • W(t) • E" 1 = e^**. 



A. Longitudinal relaxation of spins 

Now the spin excitation p(t) can be written with the 
help of Eqs. © and (|4ip in the vectorized form 



p(i) = (E-W(t)-E- 1 -C(0)) s -(h.c.) s . 



(42) 



where the subscript S means projection onto the spin 
subspace. Eq. ([42]) can be used for a fast computation. 
For the incoherent initial condition, Eqs. (fT3")l and J14j) . 
one obtains 

P(t) = ^Tr s [(E • W(t) • • (h.c.) 5 ] , (43) 

where the trace is taken over spin indices only. A more 
explicit form of Eq. (|42]) is 

N s 

p(t) - ^e-^+^cosK^-n^tlg^iU 



Ns 



Ns 



£iwC„,(o) £ L; n „c:„(o). 



(44) 



In the absence of phonon damping, T M = 0, there is a 
nonzero asymptotic value of p that can be obtained from 
the equation above by dropping all oscillating terms, i.e., 
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setting /i = v. This corresponds to the diagonal density 
matrix of the spin-phonon system. Taking into account 
L* = in the Hermitean case, one obtains 



N s 

p °° = EE i-^i 

fi n—l 



N s 

E 

n' = l 



(45) 



For the incoherent initial condition this simplifies to 



(46) 



Note that in the macroscopic limit N — > oo for a single 
spin Eq. (|44[) should assume the simple exponential form 
of Eq. ([22]) and Poo — 0. The same should be the case for 
any finite N$- 

If the phonon damping T p h is finite but small, the 
process of spin relaxation is two-stage. First, the spin 
subsystem equilibrates with the subsystem of resonant 
phonons and p(t) is mainly changing due to the time de- 
pendence of the terms with /i =^ v in Eq. (|44p , whereas 
the role of is insignificant. At the end of this stage the 
terms with /i ^ v die out, and the further slow relaxation 
is governed by T^. In particular, for the incoherent initial 
condition Eq. (|44|) at the second stage of the relaxation 
becomes 



JV S N S 
3 H n=l n' = l 



(47) 



Since for macroscopic systems the number of different 
values of in this expression is very large, the depen- 
dence p(t) is a combination of many different exponen- 
tials, i. e., p{t) is nonexponcntial. 

In the limit we study in the paper, fco r o ^ F the so- 
lution p(t) is actually the same for coherent and inco- 
herent initial conditions. The only difference is that for 
the incoherent initial condition averaging over the initial 
phases of spins leads to reproducible results for different 
realizations of the nondiagional elements of <&. For co- 
herent initial conditions, one obtains somewhat different 
results for different realizations of the spin-phonon ma- 
trix elements, related to location of the individual spins 
in space. These differences persist with increasing the 
number of spins and phonon modes, so that there is no 
self-avereaging. The computation of p(t) in the coherent 
case is faster but averaging over spin configurations is 
needed. 



that satisfies the normalization condition |co| + 
J2i \ci\ = 1 for the wave function of Eq. ([5]) in the case 
of initial phonon vacuum. In Eq. (|48]) 9 is the angle be- 
tween the spin vector and the z-axis. In the initial state 
one has 



Ns sin cos ( 



(49) 



that can be used to define the normalized transverse spin 
polarization 



(50) 



f±{t)=(<T±)t/(*±) 



that satisfies /j_(0) = 1. With the help of Eq. (gT]) one 
obtains the vectorized expression 



f±(t) 



1 

TVs 



N s 

E 

nn' — l 



E-W^-E- 1 ), 



(51) 



Note that this formula is explicitly independent of the an- 
gle 9. An alternative expression for f±(t) following from 
Eq. (|40|) has the form 



f±(t) 



5> 



where 



-l N S N S 



n=l n'=l 



In the case of undamped phonons one has L*„ 

N S 



(52) 



(53) 



Rn 



and T M = (1/N S 



is real. Then, dropping 
oscillating terms in Eq. (|5"2"|) at large times one obtains 

the asymptotic value f±(oo) = yz^j I^m| 2 [ c -f- transi- 
tion from Eq. ([44]) to Eq. (|45|) ]. It can be easily shown 



that in the presence of inhomogeneous broadening and in 
the absence of the coupling to phonons fj_(oo) — > in the 
thermodynamic limit. The same should hold in the pres- 
ence of both inhomogeneous broadening and spin-phonon 
interaction since the former should be sufficient to cause 
complete decoherence. Analysis of the principally impor- 
tant case without the inhomogeneous broadening should 
be postponed until obtaining numerical results. For a 
single spin, coupling to phonons causes decoherence with 
the rate T/2, as was stressed at the end of Sec. (jlll Al) . 
In the case of the phonon bottleneck this decoherence 
should be slowed down since a few phonon modes couple 
to many spins, still the expected result is f±(oo) = 0. 



B. Transverse relaxation of spins 

Let us now consider the time evolution of the trans- 
verse spin polarization given by Eq. (fTTj) starting from 
the fully coherent initial condition 

c 4 = —= , c = cos 9 (48) 
VNs 



C. Energy distribution of emitted phonons 

The method formulated above can be used to find the 
state of the phonon system resulting from spin-phonon 
relaxation. Obviously it can be done for undamped 
phonons only. In the case of many spins there are no 
analytical results for the occupation numbers of emitted 
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phonons rik- in the literature. However, one can express 
Uk through the eigenstates of the dynamical matrix 3? 
defined by Eqs. and ([55]) . For the wave function 
given by Eq. ([5]) one has n k = |Qc| 2 • Labeling the phonon 
modes by the discrete index I, one can express rik = rii 
through the state vector C given by Eq. ([4*0)1 . For the 
incoherent initial spin state with the help of Eq. (fl~4"|) one 
obtains, asymptotically, 



1 n s 

m(°°) = jf^2^2\ R im 

fj, ri—1 



\Rfi,N s +l\ 2 ■ 



(54) 



For Ns — 1 this expression should reproduce Eq. ([25 



VI. ANALYSIS OF THE EIGENSTATES OF THE 
SPIN-PHONON SYSTEM 

A. Spinness and off-resonance phonon emission 

Although the formalism of the preceding section is suf- 
ficient to describe the dynamics of the spin-phonon sys- 
tem in terms of transition between different bare (unper- 
turbed) modes, it is also interesting to look closer at the 
true spin-phonon eigenstates. Throughout this section 
we consider phonons as undamped, i* n = i? Am . The 
eigenstates are superpositions of spin and phonon states, 
so the first question would be to determine the fractions 
of spin and phonon states in any eigenstate, or, as it 
can be termed, their "spinness" and "phononness" . For 
instance, the spinness of the state pi is defined by 



distorted even relatively far from the resonance. This 
means, dynamically, that the number of phonon modes 
that exchange excitation with spins is not just a fixed 
number Nr determined by the single-spin relaxation rate 
r [see Eq. ([271 ] and it increases with B. Of course, Eq. 
(|57p becomes inapplicable near the resonance. However, 
one can figure out its behavior for B ^> 1. In this case 
the spinness of the states near the resonance should ap- 
proach 1. One can reproduce this behavior by adding 
(f?/2)r 2 in the denominator of Eq. (f5T[) . This is in ac- 
cordance with Eqs. ([M]) and l[55|) that define the number 
of phonon modes on speakng terms with the spins. 

For B< 1 we will see that the spinness remains small 
everywhere including the resonance. This means that 
initially prepared states with the excitation localized on 
spins will decompose over the true eigenstates that have 
a very small fraction of spin states, so that the excitation 
will migrate completely into the phonon subsystem. To 
the contrary, for B 3> 1 the initially prepared state will 
decompose over eigenstates that have spinness close to 
1, so that the excitation mostly remains on spins that is 
the phonon bottleneck. 

Using Eqs. (f5"6")) and ([54")) allows one to obtain the prob- 
ability of emission of a phonon with a wave vector k far 
from the resonance. With u => q one obtains 



N S 

5 q i=l 



V ,2 1 W*\ 2 



N s 



n ( Wk - u> y 



(58) 



N S 

Spinness ^ = \R 

n=l 



/in 



(55) 



independently of the bottleneck parameter B. This result 
is in accord with Eq. (|25p . 



while the phononness is defined by a similar expression 
with summation over phonon indices. Note that spinness 
enters the expression for the asymptotic spin excitation 
Poo, Eq. (l4*6"|) . Obviously the sum of spinness and phonon- 
ness of any state ju is 1. Spinness summed over pi gives the 
total number of spins. Far from the resonance the eigen- 
states are mainly phonon states, so that their spinness 
is small and can be calculated perturbatively. Labeling 
these states by the wave vector k instead of pi, one has 
i? qk = ^ q k, where 5 is the Kronecker symbol, and, in the 
first order of the perturbation theory, 



R 



1 



A;: 



With this one obtains 



N - 



Spinness^ 



n s \Vk\ 2 

(cj k - uoY 



B 



2 (ct> k - ui Q ) 



(56) 



(57) 



The last expression was obtained for |Vk| 2 = \V\ 2 with 
the help of Eq. ([521 . One can see that for a large bot- 
tleneck parameter the interaction of phonons with spins 
becomes large, so that phonon modes are noticeably 



B. Spin-phonon hybridization 

To quantify the spin-phonon hybridization in the gen- 
eral case, one can introduce the average phonon detuning 
from the spins u>k — u>q in any eigenstate pi as 



{LU k - LO ) 



J2k l-R/ik| 



(59) 



where the denominator is the phononness of the eigen- 
state u introduced above. Note that this definition 
does not include the momentum carried by spins. Thus 
Eq. (|59p principally cannot completely reproduce the re- 
sults for dense magnetic systems. Plotting (u>k — ^0)^ 
vs the energy (frequency) eigenvalues f2 M reveals that 
near the resonance there are pairs of the same detun- 
ings (wk-wo) p for two different frequencies f^. In 
other words, the frequencies corresponding to the same 
(wt — c<Jo) M are split because of the spin-phonon inter- 
action. Far from the resonance the relation between 
(k>k — w o) M and the frequency corresponds to a pure 
phonon mode. 
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C. Resonance scattering of phonons 

One also can study the admixture of other phonon 
modes to the given phonon mode because of the phonon 
scattering on spins. This can be described by the disper- 
sion 




(60) 

Far from the resonance the eigenstates are almost pure 
phonons, so that 5 {wk — ^o}^ is small. At resonance 
5 {^k — ^ol^ has a maximum that corresponds to the 
resonance scattering. 

VII. NUMERICAL RESULTS AND ANALYSIS 
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FIG. 1: The bottleneck plateau in the case of inhomogeneous 
broadening with the Gaussian line shape vs Sloo- 



This section has been omitted to reduce the size. 
Please, get the full text of the paper here: 

www.lehman.edu/faculty/dgaranin/Bottleneck2.pdf 



VIII. INHOMOGENEOUS BROADENING OF 
SPIN LEVELS 

In some systems the inhomogeneous broadening de- 
fined by Eq. ([3]) is much larger than the natural spin 
linewidth T or the spin-phonon splitting A. The number 
of spins within the frequency interval duo around the 
frequency ujq is 



dN s = Nsp s {too)dLo , 



(61) 



where the spin density of states p s (u>o) satisfies 
J p s (ujQ)dujQ = 1. An example is the Gaussian line 
shape 



Ps(wo) 



1 



'2ttSojq 



exp 



(u) Q - LOq) 
2 (5uj ) 2 



(62) 



In the presence of a large inhomogeneous broadening the 
number of phonons on speaking terms with spins can be 
estimated as Ns Uo ~ Nirp ph (u>o)Sujo that should replace 
Nr of Eq. (|2"T)) in the definition of the bottleneck pa- 
rameter B that becomes smaller than for the pure model 
by a factor of Sluo/T 3> 1. Consideration of this kind 
can be found in Van Vleck's original paper— and subse- 
quent publications. In fact, as we shall see below in this 
section, the relation between the numbers of spins and 
phonon modes that can exchange excitation is different 
in different frequency regions within the inhomogeneous 
spin line width Suq and a frequency-resolved description 
of the phonon bottleneck is possible. 

In simulations the macroscopic limit is achieved if the 
average distance between the neigboring spin levels be- 
comes smaller than the natural line width T 



1 



N s p s (uJo) 



«r, 



(63) 



c.f. Eq. (f2T)|) . Since both of these conditions should be 
satisfied, both N and Ns should be large enough. The 
results for the bottleneck plateau p^ vs the inhomo- 
geneous broadening Sluq in Eq. (|62[) . obtained numeri- 
cally from Eq. (|4(i>[) are shown in Fig. [TJ Since the num- 
ber of phonons that can exchange energy with spins in- 
creases with Sujo, the bottleneck plateau Poo decreases. 
For B < 1, the most pronounced decrease occurs at 
Sluq ~ r, the crossover from the natural line width to 
the inhomogeneous line width. In the case B > 1 this 
crossover occurs at Sloq ~ >/BT ~ A, in accordance with 
the comments after Eq. (|57|) . 

Clearly for Slo q 3> T, A much more phonons can 
exchange their energy with spins, so that the bottle- 
neck condition is alleviated and spin relaxation is facili- 
tated. In this case the problem simplifies since spins and 
phonons exchange energy only within the frequency in- 
terval of order max(r, A) that is much narrower than the 
inhomogeneous line width Sloq. Thus one can split the 
latter into the frequency intervals Awo around too that 
satisfy max(r, A) -C Aiv -C Sujq and consider the energy 
exchange between spins and phonons in each of these in- 
tervals independently. The spin and phonon densities 
of states within each interval Aw can be considered as 
constants and they define the bottleneck parameter in a 
frequency interval around luq 



_ N s p s (uj ) 
B«j = 



n s - 



N P P h (<*>o) P P h (wo) 



(64) 



Since the phonon density of states is a smooth function, 
one can replace p (ojq) =>■ p (ujq) . On the other hand, 
Pg(wo) and thus B U]Q have a maximum at ojq = lJq. One 
can parametrize 



B, 



Ps(wo) 



Bu, , 



(65) 



where is the bottleneck parameter at the center of 
the line. For a Gaussian line shape of Eq. (f6"2"|) one has 



B- 



ns 



27rp P h ( w o) Su) 



(66) 
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The spin excitation reaches a frequency-dependent 
plateau Poo(wn) that depends on B Wo . The average over 
the spin line shape Psi^o) of Eq. (foTj) now becomes 



P- 



daj p s {u>o)p 00 (w ). 



(67) 



Evidently poZ < poo (jjJo) since the bottleneck effect weak- 
ens away from the center of the spin band. The numeri- 
cally found dependence of Poo(wn) * s ver y close to 



(68) 



that can be expected from general statistical arguments. 
This important formula will be used below to obtain re- 
sults for the bottleneck plateau taking into account the 
inhomogeneously broadened spin line shape. The illus- 
trations will be done for the Gaussian line shape of Eq. 
161). 



In the case ifer <C 1 one obtains 



= B- 



duj{) 



As(^o) 



that with the help of Eq. (|62j) yields 



= B r 



If £w 3> 1 , one can use 



Poo = 1 - 



du>< 



PsM 

' 1 + B,„ n 



(69) 



(70) 



(71) 



that follows from Eq. (|68|) . The integrand of this expres- 
sion is close to p s (uJo)/ Bu^ = const for B UI() > 1 and it 
decays abruptly further from the center of the spin line 
where B UJa < 1. Let us define oj* that satisfies B u * — 1, 
that is, Ps(u*) = p s (uJo)/B<zq. Then from Eq. (JTT|) one 
obtains, approximately, 



^1 ^5u*p s {uj a ) 

Poo — 1 R : 



(72) 



where 6u>* 



I ■ For the Gaussian line shape of 



Eq. (f62|) one has Sui* — Sujq \/2\n B-^ and, finally, 



P- 



1 - 



In Bu 



(73) 



Fig. d shows Poo in the whole interval of B-^ calculated 
numerically from Eqs. (J62j> , ([67]), and ([65| . 

Next we calculated the time dependence of the spin 
excitation p{t) for spins within a frequency interval Ao>o 
around luq that satisfies max(r,A) <C Aluq <C Sloq, as 
explained above. Again, the incoherent initial condi- 
tion for spins leading to Eq. (|4"6")) was used. Only spins 
and phonon modes within the interval Aluq were taken 
into account while all other spins and phonon modes 
have been ignored. This allowed to greatly reduce the 
computation time. The results have been shown to be 
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FIG. 2: The bottleneck plateau in the case of strong (<5wo S> 
F, A) inhomogeneous broadening with the Gaussian line 
shape vs the bottleneck parameter B„ . 
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FIG. 3: Time dependence of the spin excitation p(t) for the 
spins with frequencies around ujq in the case of strong inho- 
mogeneous broadening. 



practically independent of Aujq as soon as the condi- 
tion r <g; Awo is fullfilled. The distribution of spin fre- 
quencies within Auo was taken equidistant (similarly to 
the phonon modes) that allowed to eleminate statistical 
scattering. In this realization of the model, spins and 
phonon modes form two equivalent groups interacting 
with each other. The results of computations for un- 
damped phonons are shown in Fig. [3J Note that here 
oscillations visible in Fig. ?? are completely washed out. 
The asymptotic values of p are in accord with Eq. (|68[) . 
Having the results for p(t) for any ajn, one could perform 
now integration over spin frequencies luq similarly to Eq. 
(|67)l using, e.g., Eq. ([6^|) . 

One can compare the results of the present calcula- 
tion within the frequency interval Aujq with the results 
of the earlier general calculation shown in Fig. [1] The 
connection is provided by the identity 



2 Sojr 



-B 



(74) 
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where p is the mass density and Vt is the speed of the 
transverse phonons. The most recent derivation of this 
formula for the transitions between any adjacent spin lev- 
els can be found in Appendix A of Ref. [lg. In zero field 
one has huj = E_ s +i - Es = (2<S - l)D and Eq. (JTSJ) 
can be cast into the elegant form 



r 



_s_<4 
i2tt nf 



m 
h 



1/4 



(76) 



The phonon density of states of Eq. (|20|) . multiplied by 
the number 2 of transverse phonon modes, has the form 



t \ - 1 

M-o)-^, 



1/3 ' 



(77) 



FIG. 4: The same in the case of damped phonons. Again, for 
B>1 the effective spin relaxation rate is much smaller than 
the phonon relaxation rate T p h. 



where vq is the unit-cell volume. The frequency fio is 
related to the Debye frequency f2_D as = {Q^Y^&d- 
Thus the bottleneck parameter B of Eq. ([28)) becomes 



that follows from Eqs. (J28J), ([64}, and ([62]). For in- 
stance, the rightmost point of the curve B = 1 in Fig. 
Q] is Poo ~ 0.070 for 5uj q /T = 30. Eq. flU} yields then 
£^ ~ 0.0418. For such small B^ one can use Eq. (|TD|) 
that yields ~ 0.030. The disagreement with the value 
Poo ~ 0.070 can be explained by the fact that the num- 
ber of spins Ns = 64 in the general calculation for B = 1 
is too small to reach the asymptotic result given by Eq. 
(155)1 and statistical scattering is still substantial. For the 
rightmost point of the curve B = 3 in Fig. Q] the dis- 
agreement between the results of the two calculations is 
smaller. 

The results for spin relaxation in the case of damped 
phonons are shown in Fig. 2J One can see that for B 3> 1 
the effective spin relaxation rate is much smaller than the 
phonon relaxation rate r p h, similarly to the case without 
inhomogeneous broadening, see Fig. ??. 



IX. IMPLEMENTATION FOR MOLECULAR 
MAGNETS 

Let us work out the general expressions and estimate 
the parameters that govern the bottlenecked spin relax- 
ation for molecular magnets, in particular, for the most 
popular compound Mni2. Historically, the phonon botte- 
neck was first observed in other systems. However, there 
is an experimental evidence of the phonon bottleneck in 
molecular magnets as well. On the other hand, molecular 
magnets are especially convenient because of the univer- 
sal form of the spin-phonon relaxation that does not de- 
pend on any unknown spin-phonon coupling constants. 15 
The spin relaxation between the adjacent levels of the 
uniaxial spin Hamiltonian —DS 2 is due to the rotation of 
the crystallographic easy axis by the transverse phonons. 
The rate of decay from the first excited state |— S + 1) to 
the ground state \—S) is given by the formula 



S{2S -l) 2 D 2 u; 3 
12nhpvf 



(75) 



B = 



12tt 2 QfCl 3 D 
~ S 



-ng. 



(78) 



For Mnia one has S = 10, p = 1.83 g/cm 3 , v = 3716 
A 3 , and from the heat-capacity measurements-ii follows 
hVL D /k B ~ 38 K thus hh D /k B ^ 10 K. Further one 
obtains v t — 2 x 10 3 m/s and Ml t /ks — 210 K, whereas 
fkoo/kB — 12 K. Plugging these parameters into Eq. ((75)) 
one obtains B ~ 5 x 10 5 ns. This means that for a non- 
diluted Mn 12 crystal, ng = 1, the bottleneck parameter 
is huge. 

Of course, for non-diluted magnetic crystals the 
physics includes the effects of coherence and it is more 
complicated than just the phonon bottleneck. The pure 
bottleneck situation is realized for a sufficient dilution, 
so that /co r o <^ 27T and the phases of emitted and re- 
absorbed phonons can be considered as random. Using 
r o = {vo/ns) 1 / 3 for the average distance between the 
neighboring magnetic molecules and ko — wo/"t, one can 
rewrite the condition of sufficient dilution as 



( 



\27rfi; 



(79) 



We call n* s critical concentration or critical dilution. 
With the parameters above, one obtains n* s ~ 0.008. 
Even at this dilution, the bottleneck parameter remains 
huge, 



B* 



i2tt 2 njn 3 D 



s 



2ttS 



(80) 



that numerically yields B* ~ 4000. Note that for spin 
transitions between excited levels the energy differences 
Hujq = E m+ i — E rn = — (2m + 1) D are smaller than 
above, thus the values of B are even larger. 

Now from Eq. (|19p one obtains the estimation of the 
spin-phonon matrix element 



W\ 2 = 



r 



s uj 3 n n 3 



27T P ph (uJ Q ) 



24 



(81) 



11 



(that also could be obtained directly!). This yields the 
spin-phonon splitting of Eq. (|31[) in the form 



parameter of Eq. With the help of Eq. |(77|) 

one obtains 



S luqVI'j^ 



and 



A* = 




(82) 



(83) 



Numerically one obtains hA/ks — 40 v / ns" mK and 
HA*/k B ~ 3.6 mK. 

Let us now discuss the role of inhomogeneous broad- 
ening on the phonon bottleneck in molecular magnets. 
First, there is the dipole-dipole interaction (DDI) that 
is as strong as about EoDi/ks — 67 mK between the 
two neighboring Mni2 molecules. This would result in 
the change of the spin transition frequency for adjacent 
spin levels by {Eddi /ks) /S ~ 6.7 mK. Sometimes, for 
simplicity, the DDI is considered as a kind of inhomoge- 
neous broadening. This is an oversimplification, at least 
for weakly excited states considered here. Rigorous treat- 
ment of the Landau-Zener effect at fast sweep with an ac- 
count of both DDI and true inhomogeneous broadening 
(random hyperfine fields}^ shows that these two effects 
compete with each other, rather than simply add. In 
the present case, the magnetostatic field smoothly varies 
along the crystal (if the crystal shape is non-elliptic), so 
that its change on the unit-cell distance is very small. 
Thus the DDI as a source of inhomogeneous broadening 
will be neglected. 

The greatest source of inhomogeneous broadening in 
Mni2 is the hyperfine interaction with their own Nj — 
12 nuclear spins I = 5/2. With the hyperfine coupling 
A between the total electronic spin S and each of the 
nuclear spins I of A/k B = 2 mK (Ref. QJ), the dispersion 
of the hyperfine field SHhf on the electronic spin is given 



SH HF = ^, ai = ^l (I+1): (84 ) 

whereas the inhomogeneous line shape is Gaussian. Nu- 
merically one obtains SHhf — 8.8 mT. For the transi- 
tion between adjacent spin levels for the inhomogeneous 
broadening Sluq in Eq. (|62j) one has HSluq = gu B SHnF, 
thus hdcoo/kB — 12 mK. In another popular compound 
Fes the inhomogeneous broadening is mainly due to the 
DDI with nuclear spins of hydrigen atoms present in mag- 
netic molecules, SH ~ 0.8 mT according to Ref. l2ll . 

One can see that in Mni2 the inhomogeneous broaden- 
ing 12 mK is smaller than the spin-phonon gap hA/ks — 
40 mK in the non-diluted case, thus it can be neglected 
in the first approximation. On the other hand, as said 
above, the pure bottleneck case requires dilution below 
n* s of Eq. (|79|) . and for ns < n s one has Sujq > A. For Fe s 
the inhomogeneous broadening begins to play a role for 
a very strong dilution where measurements are difficult. 
If the inhomogeneous broadening is dominating, the re- 
laxation is governed by the frequency-interval botleneck 



_ ns Tr 2 n 3 D 

— c 

At the critical dilution this becomes 



n% Tr 2 Q 3 n 

j^* o JJ_ 

U ° V27T UJqSlUo 



27/2^3/2^ - U - U1D ^o' 



(85) 



(86) 



That is, for any line with a well-defined central frequency 
uJo the bottleneck parameter is large. With the above 



parameters for Mni2 one obtains Bjz 
U^/Slu ~ 1000, and B± 



: 2 x 10 d n s , 
~ 16 that is still large. Then 
Eq. d73j) yields 1 - ~ 0.11, that is, only a 10% of 
the initial spin excitation can relax before the bottleneck 
plateau is reached. Thus the inhomogeneous broadening 
in Mn 12 does not fully resolve the PB, and in Fe 8 its 
effect is even smaller. 



X. DISCUSSION 

The problem of the phonon bottleneck considered in 
this paper in the weak-excitation limit is only a part of a 
larger problem of collective spin-phonon relaxation. The 
bottleneck in the pure form occurs for magnetically di- 
luted systems that satisfy the condition k^r^ 1, where 
ko is the wave vector of a resonant phonon and ro is the 
typical distance between the neighboring spins. If this 
condition is violated, one cannot consider the phases of 
emitted phonons reaching other spins as random, and the 
interference effects become important. The well-known 
example of interference effects is superradianc o 3 i 12 that 
requires, among other conditions, a coherent initial con- 
dition of spins. In contrast to the bottleneck, superradi- 
ance dramatically increases the relaxation rate, so that 
the two effects should compete. On the other hand, de- 
structive interference effects in the case of incoherent ini- 
tial condition or dynamical loss of coherence due to the 
inhomogeneous broadening can lead to suppression of re- 
laxation that resembles the PB but has a different physi- 
cal origin. For this reason, using the results of this paper 
to interpret experiments should be done with care as the 
experimental observations can be a mixture of different 
effects. In particular, the results obtained in this paper 
are not applicable to non-diluted molecular magnets. 

Considering the weak-excitation limit in this paper al- 
lowed to drastically simplify the Schrodinger equation 
and obtain numerically exact results for the PB effect 
with and without inhomogeneous broadening of spins and 
phonon damping. It was confirmed that the bottleneck 
parameter B quantifying the statistical weights of spins 
and resonant phonons, introduced in Ref. 11, plays the 
main role in the problem. Fundamentally the most in- 
teresting case corresponds to the pure model without in- 
homogeneous spin broadening and phonon damping. For 
this model, similarly to the results of Ref. [ll], the spin 
excitation p(t) was shown to oscillate approaching the 
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bottleneck plateau, see Fig. ??. However, these oscilla- 
tions have a smaller amplitude and are stronger damped 
than the analytical results of Ref. [ll|. The frequency of 
these oscillations corresponds to the gap A between the 
two branches of the hybridyzed magnetoelastic waves at 
resonance [see Eq. (j3"U|) and Fig. ??] 

Here it was shown that the bottleneck parameter effec- 
tively decreases in the presence of inhomogeneous broad- 
ening that alleviates the bottleneck condition. If the in- 
homogeneus spin line width Sojq exceeds the spin-phonon 
gap frequency A/h, the splitting of spin and phonon 
modes is not resolved and the oscillations of p(t) are 
washed out, see Fig. [3] 

Inclusion of the ad hoc phonon relaxation rate F p h in 
the theory describes the second, post-plateau, stage of 
the spin relaxation. An important observation is that 
for B 1 the corresponding relaxation rate is much 
smaller that T p h (see Figs. ?? and rjj and the spin re- 
laxation is non-exponential [see Eq. ([4"?)) and Fig. ??]. 
One should stress, however, that it is not completely sat- 
isfactory to plug an ad hoc phonon relaxation into the 
theory. The latter should be treated in this case as taken 
from the experiment. This can lead to a problem since 
one of the main sources of the observed phonon damp- 
ing can be their scattering on spins that is already taken 
into account by the very spin phonon interaction, Eq. 
(TJJ. On the other hand, this and other kinds of elastic 
scattering cannot help the system to reach complete equi- 
librium since these processes conserve the energy and do 
not transfer excitation from the narrow group of resonant 
phonons to the rest of the phonon bath. 



It would be of a principal importance to generalize the 
theory of phonon bottleneck for the higly excited initial 
states of the spin subsystem. However, the Schrodinger 
equation in this case is not amenable to a direct numerical 
solution. In Ref. [l^ it was argued that, since collective 
motion of spins in the regions large in comparizon to the 
typical distance between the neighboring spins involves a 
large number of atomic spins, the problem can be consid- 
ered classically. Thus the spin and phonon operators in 
Eqs. (9) and (10) of Ref. Qjl had been replaced by classi- 
cal variables. This gives a possibility to numerically treat 
highly excited states of the spin-phonon system without 
the combinatorial explosion of the full SE. On the other 
hand, linearization of the classical equations of motion 
of Ref. [H near the ground state yields equations that 
are equivalent to the truncated SE, Eq. ([7]), in the low- 
excitation limit. This suggests that classical equations 
Ref. [l2] are indeed a good approximation for the spin- 
phonon problem in the whole energy range. It would be 
interesting to investigate what are the quantum correc- 
tions to these equations. 
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